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Summary 


The  following  objectives  have  been  achieved: 

•  developed  a  method  for  producing  molecular  oxygen  in  its  first  electronically  excited  state 
(designated  CMa'AgJ),  based  on  the  technique  pioneered  by  Dr  A.  A.  Viggiano  at  AFRL 

•  produced  an  absolute  calibration  of  the  02(a)  concentration  using  1270  nm  emission  from 
02(a*Ag  -  the  yield  of  02(a)  from  CD  in  the  generator  ranged  from  16  to  26% 

•  studied  the  reaction  kinetics  at  300  K  of  atomic  Mg,  Ca  and  Fe  with  02(a)  in  a  fast  flow  tube, 
where  the  metal  atoms  were  produced  either  by  thermal  evaporation  (Mg  and  Ca)  or  pulsed 
laser  ablation  (Fe),  and  then  detected  by  laser  induced  fluorescence 

•  carried  out  high  level  electronic  structure  calculations  to  explore  the  potential  energy  surfaces 
for  these  metal  atom  +  02(a)  reactions,  as  well  as  the  reaction  SiO  +  02(a) 


Conclusions: 

•  the  reaction  Ca  +  02(a)  mostly  occurs  via  recombination  to  produce  Ca02(*Ai),  with  a  rate 
coefficient  that  is  ~80  times  larger  than  for  the  reaction  Ca  +  02(X).  There  is  also  limited 
switching  through  a  conical  intersection  between  singlet  and  triplet  surfaces,  leading  to  the 
bimolecular  products  CaO  +  O 

•  the  reaction  Mg  +  02(a)  reaction  occurs  exclusively  by  recombination  on  a  singlet  surface, 
producing  Mg02(*Ai).  The  rate  coefficient  is  -4000  times  larger  than  for  the  reaction  Mg  + 
02(X).  However,  there  are  no  surface  crossings  between  the  initial  singlet  and  lower-lying 
triplet  surfaces 

•  the  reaction  Fe  +  02(a)  produces  FeO  +  O,  although  with  a  probability  of  only  -  0.1%.  There 
is  no  evidence  for  recombination,  suggesting  that  this  reaction  proceeds  mostly  by  near¬ 
resonant  electronic  energy  transfer,  producing  02(X)  and  Fe(a5F)  (the  lowest-lying  excited 
state  of  Fe) 

•  the  reaction  SiO  +  02(a)  has  a  significant  barrier  of  55  kJ  mol1  and  also  requires  a  singlet- 
triplet  surface  crossing  in  order  to  produce  Si02  +  O 

•  the  reaction  Si+  +  02(a)  (measured  by  Dr  Viggiano  at  AFRL)  is  the  most  important  removal 
process  for  Si+  ions  in  the  daytime  atmosphere  between  85  and  107  km 

•  the  reaction  of  Ca  with  02(a)  is  a  factor  of  -4  times  slower  than  the  reaction  of  Ca  with  O3, 
during  daytime  between  85  and  110  km 

•  the  reactions  of  Fe,  Si,  SiO  and  Mg  with  02(a)  are  too  slow  to  be  atmospherically  significant. 
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kinetics  of  Ca,  Mg  and  Fe  with  02(a  1  AgJ,  J.  Phys.  Chem.  A,  to  be  submitted  spring  2012. 
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1.  Introduction 


The  first  excited  state  of  molecular  oxygen,  (Via1  A„),  is  produced  in  the  mesosphere-lower 
thermosphere  (MLT)  region  by  photolysis  of  O3  at  wavelengths  shorter  than  320  nm: 

03  +  hv  O('D)  +  OoCa'Ag) 

(Tfa1  Ag)  is  comparatively  long-lived  in  the  MLT.  The  quenching  life-time  is  more  than  4  hours, 
much  longer  than  the  phosphorescence  lifetime  of  72  minutes  for  emission  at  1.27  pm  in  the  infra¬ 
red.  Since  the  rate  of  O3  photolysis  in  the  MLT  is  ~8  x  10  3  s'1,  then  during  daytime  the  steady-state 
ratio  [  02(a 1  Ag) ]/[ O3J  is  about  30  [Brasseur  and  Solomon,  2005]  and  the  daytime  concentration  of 
CLCa'Ag)  around  90  -  100  km  is  ~  5  x  109  cm  3  [Batista  et  al.,  1996].  After  sunset  CL^Ag)  decays  by 
an  order  of  magnitude  every  2.8  hours.  CL^Ag)  contains  almost  1  eV  of  electronic  excitation 
compared  with  ground-state  02(XaZg"),  and  thus  has  the  potential  to  be  significantly  more  reactive, 
although  it  is  not  a  radical  species. 

The  purpose  of  this  project  was  to  explore  the  impact  of  CL^Ag)  on  the  constituents  of  the  MLT 
which  are  produced  by  meteoric  ablation.  The  interactions  of  metallic  and  silicon  species  with 
02(a1Ag)  have  not  before  been  studied  (to  the  best  of  our  knowledge).  However,  following  the  recent 
development  of  a  new  method  for  producing  CLl  'Ag)  cleanly  and  at  significant  concentrations  [Midey 
et  al.,  2008],  such  a  study  has  become  feasible.  The  work  was  undertaken  in  collaboration  with  the 
research  group  of  Dr  A.  A.  Viggiano  (Air  Force  Research  Laboratory,  Space  Vehicles  Directorate), 
who  pioneered  this  new  technique  of  producing  CL^Ag)  for  laboratory  kinetic  studies  [Midey  et  al., 
2009], 

The  major  source  of  metals  in  the  earth’s 
upper  atmosphere  is  the  ablation  of  the 
roughly  50  tonnes  of  interplanetary  dust  that 
enters  each  day  from  space  [Plane,  2003]. 
This  gives  rise  to  the  permanent  layers  of 
metal  atoms  that  occur  globally  in  the  MLT 
region  between  about  75  and  110  km 
(Figure  1).  Four  metals  -  Na,  Fe,  K  and  Ca 
-  have  been  studied  intensively  during  the 
last  two  decades  using  the  ground-based 
resonance  lidar  technique  [Bowman  et  al., 
1969;  Eska  et  al.,  1999;  Gerding  et  al., 
2000;  Kane  and  Gardner,  1993].  Metallic 
(and  silicon)  ions  have  been  detected  by 
rocket-bome  mass  spectrometry 
[Grebowsky  and  Aikin,  2002;  Kopp,  1997]. 
Most  recently,  satellite -borne  spectrometric 
observations  of  the  earth’s  dayglow  have 
been  used  to  retrieve  the  global  layers  of 
metal  atoms  and  ions,  in  particular  Mg  and 
Mg+  [Fan  et  al.,  2007;  Scharringhausen  et 
al.,  2008]. 

Meteoric  ablation  injects  metal  atoms  and  ions  directly  into  the  MLT  [Vondrak  et  al.,  2008].  These 
species  then  undergo  neutral  and  ion-molecule  chemistry.  Below  90  km,  metal  atoms  are  rapidly 
oxidised  by  a  series  of  reactions  involving  O3,  CL,  FLO  and  CO2  to  reservoir  species  such  as  oxides, 
hydroxides  and  bicarbonates.  However,  these  species  react  with  either  atomic  O  or  H  (and,  to  a  lesser 
extent,  undergo  daytime  photolysis)  to  regenerate  metal  atoms.  Thus,  the  atomic  metal  layers  only 
appear  above  80  km,  where  atomic  O  and  H  increase  by  several  orders  of  magnitude  and  are  present 


Figure  1.  Vertical  profiles  of  the  annual  mean 
concentrations  at  northern  hemisphere  mid-latitudes  of: 
Na  (Plane,  J.  M.  C.  et  al.  J.  Geophys.  Res.  104,  3773- 
3788,  1999);  Fe  (Helmer,  M.  et  al.  J.  Geophys.  Res.  103, 
10913-10925,  1998.);  K  (Eska,  V.  et  al.  J.  Geophys.  Res. 
104,  17173-17186,  1999);  and  Ca  (Gerding,  M.  et  al.  J. 
Geophys.  Res.  105,  27131-27146,  2000). 
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both  during  day  and  night  [Plane,  2003].  This  ensures  rapid  recycling  of  the  metal  atoms  from  their 
reservoirs.  The  small  scale  height  of  2  -  3  km  (the  scale  height  is  the  distance  over  which  the 
concentration  changes  by  a  factor  of  e)  on  the  undersides  of  the  layers  (Figure  1)  closely  follows  the 
atomic  O  fall-off.  Below  about  85  km  the  reservoir  species  are  permanently  lost  by  condensation 
onto  meteoric  smoke  particles.  These  are  approximately  1  nm-size  particles  which  form  from  the 
polymerization  of  metal  oxides,  carbonates  and  silicates  produced  by  meteoric  ablation  [ Bardeen  et 
al.,  2008;  Hunten  et  al.,  1980;  Kalashnikova  et  al.,  2000;  Saunders  and  Plane,  2006].  Meteoric 
smoke  particles  also  probably  provide  ice  nuclei  for  the  formation  of  noctilucent  clouds  (NLCs)  in 
the  summer  high-latitude  mesosphere  [ Bardeen  et  al.,  2008;  Megner  et  al.,  2006].  The  ice  particles  in 
these  clouds,  with  radii  of  -50  nm,  remove  metallic  species  even  more  efficiently  [ Gardner  et  al., 
2005;  Plane  et  al.,  2004].  Smaller  ice  particles  (5  -  10  nm  radius)  affect  the  plasma  locally  and  give 
rise  to  very  strong  radar  echoes  known  as  polar  mesospheric  summer  echoes  (PMSE)  around  85  km 
[Megner  et  al.,  2006].  NLCs  appear  increasingly  to  be  a  sensitive  indicator  of  climate  change  in  the 
mesosphere  [Shettle  et  al.,  2009]. 

Above  95  km  in  the  lower  E  region  of  the  ionosphere,  there  is  a  marked  increase  in  the  concentration 
of  plasma,  which  consists  mostly  of  NO+,  C>2+,  and  electrons  [Plane,  2003].  Because  of  their  low 
ionization  potentials,  metal  atoms  charge  exchange  readily  with  these  positive  ions  to  produce 
metallic  ions.  Metallic  ions  then  form  clusters  with  No,  O2,  CO2  and  FLO  to  form  strongly-bound 
cluster  ions  which  undergo  dissociative  recombination  with  electrons  to  yield  neutral  metal  atoms. 
This  process  is  interrupted  by  atomic  O,  which  tends  to  reduce  cluster  ions  back  to  metal  ions  [Plane, 
2003].  Fe+,  Mg+  and  Si+  are  the  major  components  of  sporadic  E  layers  [  G  re  how  sky  and  Aikin, 

2002].  These  are  thin  layers  of  concentrated  plasma  (essentially  metallic  ions  and  electrons)  that 
occur  in  the  lower  thermosphere  between  about  90  and  120  km.  They  have  a  significant  impact  on 
radio  communications,  both  by  facilitating  over-the-horizon  HF  communication  and  by  obscuring 
space-to-ground  communications.  The  lifetime  of  a  sporadic  E  layer  is  largely  controlled  by  ion- 
molecule  chemistry,  and  can  vary  from  several  days  at  105  km  to  only  a  few  minutes  at  90  km 
[Woodcock  et  al.,  2006]. 


2.  Specific  objectives  of  the  project 


The  project  involved  the  kinetic  study  of  a  set  of  reactions  between  various  meteoric  species  and 
02(a).  These  reactions,  which  are  listed  in  Table  1,  are  all  exothermic  (or  very  close  to 
thermoneutral),  and  could  thus  be  rapid  and  play  important  roles  in  the  MLT.  The  reactions  are  also 
of  fundamental  interest  for  understanding  the  role  of  electronic  excitation  and  spin  in  controlling 
chemical  reactions. 

The  reaction  between  Ca  and  02(a)  (reaction  1)  could  provide  an  alternative  daytime  route  to 
oxidizing  Ca  via  O3  [Helmer  et  al.,  1993],  although  reaction  1  is  formally  spin-forbidden. 
Alternatively,  recombination  to  form  the  most  stable  isomer  of  Ca02  is  spin-allowed,  and  may  be 
considerably  faster  than  the  recombination  of  Ca  and  ground-state  02(X),  which  has  a  barrier  of 
about  6  kJ  mol1  [Campbell  and  Plane,  2001]. 

Although  not  specified  in  the  original  proposal,  we  added  a  study  of  the  reaction  between  Mg  and 
02(a)  (reaction  2).  This  reaction  contrasts  with  reaction  1  because  if  reaction  occurs  then  it  must 
involve  recombination  to  Mg02  (this  is  because  the  Mg-0  bond  is  relatively  weak  and  so  formation 
of  MgO  +  O  is  endothermic  by  164  kJ  mol1  [hide,  2006]).  Note  that  both  the  Ca  and  Mg  reactions 
could  also  quench  the  02(a)  to  the  ground  state  02(X),  although  this  would  involve  crossing  from  the 
reactant  singlet  surface  onto  the  product  triplet  surface. 

Fe  is  oxidized  to  FeO  at  every  collision  with  O3  [Helmer  and  Plane,  1994b].  Given  that  the  [02(a)]  / 
[O3]  ratio  in  the  mesosphere  is  -30,  the  reaction  between  Fe  and  02(a)  (reaction  3)  could  therefore 
substantially  increase  the  daytime  oxidation  rate  of  Fe.  This  reaction  could  also  involve 
recombination  on  a  quintet  surface  to  yield  Fe02  (or  the  more  stable  OFeO  isomer).  In  contrast,  the 
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reaction  between  Fe  and  ground-state  CF(X  Sg")  can  only  produce  FcCK  but  has  a  large  electronic 
barrier  of  about  17  kJ  mol"1  [Helmer  and  Plane ,  1994a].  Another  interesting  possibility  is  near¬ 
resonant  electronic  energy  transfer  to  yield  Fe(a5F)  +  C>2(X).  Fe(a5F)  is  the  first  electronically-excited 
state  of  Fe. 

Reaction  4  may  provide  an  important  route  in  the  mesosphere  for  oxidizing  SiO  to  SiCK  which  is  the 
most  stable  form  of  silicon  and  a  likely  precursor  of  meteoric  smoke  [Saunders  and  Plane ,  2006]. 
One  pathway  is  the  reaction  between  SiO  and  O3:  however,  we  have  shown  recently  that  this  reaction 
is  very  slow  [Gomez  Martin  et  al,  2009b],  so  that  reaction  2  could  oxidize  SiO  up  to  15,000  times 
faster,  although  it  would  involve  a  singlet-triplet  surface  crossing. 

During  the  course  of  this  project,  the  reaction  between  Si+  and  02(a)  (  reaction  5)  was  studied 
experimentally  at  AFRL  by  Dr  A  A  Viggiano’s  group  [Eyet  et  al.,  2010].  The  reaction  was  shown  to 
be  relatively  fast  (^(200  K)  =  4  x  10  11  cnr1  molecule"1  s"1).  The  atmospheric  implications  were  then 
studied  as  part  of  this  project. 

The  original  proposal  also  listed  the  reaction  between  Si  and  02(a)  for  study.  This  reaction  is 
exothermic  and  spin-allowed,  producing  SiO  +  O('D)  or  SiO(3n)  +  0(3P)).  However,  we  have 
shown  recently  that  the  reaction  between  Si  and  ground-state  02(X)  is  fast  ((k(200  K)  =  1.3  x  10"10 
cm3  molecule'1  s"1)  [Gomez  Martin  et  al.,  2009a],  so  that  the  reaction  with  02(a)  will  have  an 
insignificant  effect  on  the  oxidation  of  Si  in  the  MLT.  This  reaction  was  therefore  replaced  in  the 
project  with  a  study  of  reaction  2. 


Table  1.  Reactions  of  meteoric  species  with  02(a1Ag)  that  could  play  important  roles  in  the 
mesosphere  and  lower  thermosphere 


No. 

Reaction 

A//0  / 

kJ  mol"1 

Electron  spin 
conserved? 

1 

Ca  +  02(a1Ag)  -+  CaOl1!)  +  0(3P) 

0a 

No 

->  Ca02(1A1) 

-279  b 

Yes 

-+  Ca  +  02(32V) 

-93  a 

No 

2 

Mg+  02(a‘Ag)  -+  Mg02(1Ai) 

-149  b 

Yes 

^Mg  +  02(32:g") 

-93  a 

No 

3 

Fe  +  CMa1  Ag)  -+  FeO(5Z)  +  0(3P) 

-2+8 

Yes 

->  FeCM'AO 

-376  c 

Yes 

— >  Fe(a5F2)  +  O2 

0a 

Yes 

-+  Fe(a5D)  +  02 

-94  a 

Yes 

4 

SiO  +  02(a1Ag)  -+  Si02  +  0(3P) 

-31  b 

No 

->  OS^AO 

-273  b 

Yes 

5 

Si+  +  02(a‘Ag)  -+  SiO+  +  O 

-67  b 

Yes 

a  Using  experimental  energies  [Lide,  2006].  b  Calculated  at  the  CBS-Q  level  of  theory  [Frisch  et  al., 
2009]; c  Calculated  at  the  B3LYP/6-311+g(2d,p)  level  of  theory  [Frisch  et  al.,  2009] 
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3.  Experimental  technique 


Figure  2  is  a  schematic  diagram  of  the  fast  flow  tube  apparatus  that  was  used  to  study  the  kinetics  of 
the  Ca,  Mg  and  Fe  reactions.  The  stainless  steel  flow  tube  has  an  internal  diameter  of  37.5  mm  and 
consists  of  sections  of  tube,  cross-pieces  and  nipple  sections  connected  by  conflat  flanges  sealed  with 
copper  gaskets.  The  tube  has  a  total  length  of  1 130  mm  from  the  upstream  entry  point  of  the  carrier 
gas  to  the  downstream  laser  induced  fluorescence  (LIF)  detection  cell.  Calcium  atoms  were 
produced  continuously  by  heating  calcium  pellets  (Aldrich,  99  %)  to  1070  -  1120  K.  Magnesium 
atoms  were  produced  by  heating  magnesium  pellets  (Aldrich,  99.5  %)  to  a  temperature  between  700 
and  800  K.  The  pellets  were  located  in  an  aluminium  oxide  crucible  placed  inside  a  tungsten  basket 
heater,  positioned  1120  mm  upstream  of  the  LIF  cell.  The  Ca  or  Mg  atoms  were  entrained  in  the 
main  carrier  flow  of  No,  which  entered  the  tube  upstream  of  the  crucible.  Ca  was  detected  by  LIF  at 
422.7  nm  (Ca(41Pi  -  4IS))  using  a  Nd-YAG-pumped  dye  laser  (pulse  rate  10  Hz;  pulse  energy  10 
mJ).  Mg  was  detected  by  LIF  at  285.2  nm  (Mg(31Pi  -  3 1  So)),  after  frequency-doubling  the  dye  laser 
using  a  BBO  crystal. 


Figure  2.  Schematic  diagram  of  the  02(a)  generator  coupled  to  a  fast  flow  with  laser  induced 
fluorescence  detection  for  studying  metal  atom  reactions  (exemplified  by  Mg)  with  CLla'Ag) 

The  source  of  Fe  atoms  in  the  flow  tube  was  the  pulsed  ablation  of  a  pure  iron  rod,  using  a  Nd:YAG 
laser  (k  =  532nm,  pulse  energy  =  22-31  mJ,  repetition  rate  =  8  Hz).  The  rod  was  coupled  to  a 
stepper  motor  (via  a  vacuum  feedthrough  in  a  side-arm  of  the  flow  tube),  so  that  the  rod  could  be 
rotated  (2-4  Hz)  and  also  translated  slowly.  This  ensured  that  a  fresh  surface  of  the  rod  was 
presented  to  each  successive  laser  shot,  in  order  to  keep  the  resulting  pulses  of  ablated  Fe  as  uniform 
as  possible.  The  iron  rod  was  long  enough  (~  5  cm)  to  project  across  the  central  axis  of  the  tube. 

The  laser  was  loosely  focused  onto  the  rod  through  an  orthogonal  side-arm,  so  that  the  point  of 
ablation  was  in  the  centre  of  the  flow  tube.  The  pulse  of  Fe  atoms  was  then  entrained  in  the  No  carrier 
gas  and  transported  downstream  to  the  LIF  cell,  where  the  Fe  was  detected  by  LIF  at  248.3  nm 
Fe(x5F5  a5D4). 
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The  reactant  flow  (Oo(a)  in  He)  was  injected  via  a  side  port  downstream  of  the  crucible  (Figure  2). 
The  gas  flow  exited  the  tube  through  a  throttle  valve  to  a  booster  pump  backed  by  a  rotary  pump, 
providing  a  volume  displacement  rate  of  1 10  1  s'1.  Typically,  a  total  gas  flow  rate  of  3200  seem  was 
used  with  pressures  ranging  from  1  to  10  Torr.  The  Reynolds  number  was  always  below  80,  ensuring 
laminar  flow  within  the  tube. 

02(a1Ag)  was  prepared  using  the  technique  pioneered  by  Viggiano’s  group  [Midey  et  al.,  2008;  Midey 
et  al.,  2007;  Midey  et  al.,  2009],  where  Cl2  is  bubbled  through  a  chilled  alkaline  solution  of  H2O2: 

Cl2  (g)  +  H2O2  (aq)  +  2KOH(aq)  ^  02(a'  Ag)  /  02(X3Eg')  (g)  +  2KC1  (aq)  +  H20  (6) 

This  produces  02(a1Ag)  at  up  to  30%  yield  [Midey  et  al.,  2008].  The  02(a)  generator  constructed  for 
this  project  is  illustrated  schematically  in  Figure  2.  The  two  traps  and  the  cell  used  to  monitor 
emission  at  1270  nm  from  02(a)  were  contructed  from  Pyrex  glass.  The  first  trap  is  where  reaction  6 
takes  place.  This  trap,  containing  60  ml  of  35%  w/w  H2O2,  was  held  at  -21°C  (provided  by  a  water 
ice,  dry  ice  and  NaCl  slush).  Using  a  burette,  40  ml  of  chilled  4.04  M  KOH  was  added,  and  the  trap 
was  then  attached  to  the  rest  of  the  generator.  A  10%  CF/He  was  then  bubbled  through  the  resulting 
slush,  at  flow  rates  up  to  100  seem.  The  second  trap,  held  at  -70°C  (using  a  30%  H20/CH30H  slush), 
was  used  to  dry  the  gas  flow  by  freezing  out  H20.  The  CF/Hc  flow  then  entered  the  cylindrical 
optical  emission  cell  (length  =  100  mm,  radius  =  10  mm),  before  flowing  through  a  Teflon  valve  into 
the  fast  flow  tube. 

The  weak  emission  at  1270  nm  from  02(a-X)  emission,  exiting  through  a  window  at  one  end  of  the 
cell,  was  focused  by  a  lens  (bi-convex  BK-7  lens,/=  30  mm))  through  an  interference  filter  (centre 
wavelength  =  1270  nm,  FWHM  =  4.2  nm,  peak  transmission  =  33%)  into  a  glass  fibre  optic  bundle 
of  length  610  mm  (Figure  2).  The  light  exiting  the  bundle  was  then  focused  by  a  second  lens  onto  an 
InGaAs  photodiode  detector  (Oriel,  Model  71671).  The  photodiode  current  was  read  by  a  pico- 
ammeter  (Keithley).  The  absolute  sensitivity  of  this  optical  assembly  -  light  collection,  transmission 
through  the  fibre  optic,  and  detection  -  was  calibrated  using  a  radiometric  calibration  standard 
(Ocean  Optics  LS-l-CAL  NIST-traceable  light  source). 


Experimental  Results 


Calibration  of  the  02(a)  flow 

Figure  3  shows  that  the  InGaAs  photodiode  current  resulting  from  1270  nm  emission  in  the  optical 
cell  was  always  proportional  to  the  Cl2  concentration  entering  the  first  trap  of  the  02(a)  generator. 
This  implies  that  a  constant  fraction  of  the  02  produced  by  reaction  6  was  in  the  a*Ag  state.  The 
calibration  of  the  absolute  02(a)  concentration  was  achieved  in  two  stages: 

1.  A  computer  ray-tracing  model  was  developed  to  calculate  the  total  collection  efficiency  of  1270 
nm  photons  emitted  in  the  gas  cell.  This  model  assumes  that  02(a)  has  a  uniform  concentration  in  the 
emission  cell,  which  should  be  the  case  given  its  long  lifetime  (rrad  =  4200  s;  quenching  in  He  is 
negligible).  The  model  then  determines  the  probability  of  a  photon  emitted  at  each  point  in  the  cell 
being  captured  by  the  bi-convex  lens  and  focused  through  the  interference  filter  onto  the  entrance  of 
the  fibre  optic  bundle.  The  total  number  of  photons  entering  the  bundle  is  then  computed  by 
integrating  over  the  cell  volume.  A  correction  is  then  applied  because  the  02(a-X)  vibrational  band 
structure  is  significantly  broader  that  the  bandpass  of  the  interference  filter:  only  20.1%  of  the  total 
emission  spectral  intensity  falls  within  the  filter  bandpass.  The  result  is  that  if  the  concentration  of 
02(a)  in  the  cell  were  1.0  x  1016  cm'3,  then  the  02(a-X)  emission  power  entering  the  fibre  bundle 
would  be  0.12  nW. 

2.  With  the  light  collecting  system  removed  from  the  optical  emission  cell,  the  photodiode  current 
was  measured  as  a  function  of  the  distance  between  the  radiometric  calibration  standard  and  the 
collection  lens.  This  showed  that  the  calibration  factor  was  103  pA  nW"1.  Hence,  a  concentration  of 
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1.0  x  1016  cm"3  02(a)  in  the  cell  would  produce  a  photodiode  current  of  12.4  pA.  That  is,  the 
calibration  factor  was  8.1  x  1014  cm  3  pA"1. 

The  right-hand  ordinate  in  Figure  3  shows  the  resulting  02(a)  concentrations,  calculated  by  applying 
this  calibration  factor  to  the  detector  current  on  the  left-hand  ordinate.  The  selection  of  experimental 
runs  shown  in  Figure  3  covers  the  range  of  conversion  efficiencies  of  CI2  into  02(a)  observed  during 
the  project,  which  ranged  from  16  to  26%.  This  variation  in  efficiency  seemed  most  likely  due  to  the 
cleanliness  of  the  Pyrex  glass  traps. 
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Figure  3.  02(a)  emission  current  measured  with  the  In-Ga-As  detector  at 
1270  nm,  as  a  function  of  [CI2]  in  the  generator.  The  corresponding 
calibrated  [02(a)]  is  shown  on  the  right-hand  ordinate.  Data  from  a 
selection  of  experimental  runs  over  several  months  shows  that  the 
efficiency  for  02(a)  production  ranged  from  16  -  26%  of  the  CI2. 


Kinetic  experiments 

Reaction  rate  coefficients  were  determined  using  a  procedure  we  have  described  in  detail  elsewhere 
[Broadley  et  al.,  2007].  Taking  reaction  1  as  an  example,  the  loss  of  Ca  by  diffusion  to  the  flow  tube 
walls  and  reaction  with  02(a),  can  be  described  by  a  first-order  decay  coefficient,  k',  since  [02(a)] 
»  [Ca]: 

kf=  Wa+^°2(a)]  (I) 


where  kdiff.Ca  describes  the  loss  of  Ca  by  diffusion  and  k  is  the  rate  coefficient  for  reaction  1  (which 
may  depend  on  pressure).  An  experiment  is  carried  out  by  varying  [02(a)]  while  keeping  the  total 
mass  flow  rate  and  pressure  in  the  flow  tube  constant.  This  means  that  kdiff.Ca  is  constant,  as  well  as 
the  reaction  time  t  between  the  point  of  injection  of  the  02(a)  and  the  downstream  LIF  cell.  Since  the 
removal  of  Ca  is  pseudo  first-order, 


In 


[Ca];, 

jcar02(a)y 
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ln[Ca]rel 
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where  [Ca]g  is  the  concentration  at  the  LIF  detection  cell  in  the  absence  of  02(a),  [Ca]()  (a)  is  the  Ca 

concentration  at  the  LIF  detection  cell  when  02(a)  is  added,  and  [Ca]rei  is  the  ratio  of  these 
concentrations.  Plots  of  ln[Ca]rei/  t  versus  [02(a)]  are  shown  in  Figure  4,  for  a  range  of  pressures  in 
the  flow  tube.  The  linear  dependence  expected  from  equation  II  is  observed,  and  the  slope  of  each 
plot  gives  the  second-order  rate  coefficient  k. 


[02(a)]  / 1 012  molecule  cm"3 

Figure  4.  Kinetic  plots  showing  the  first-order  removal  rate  of  Ca  as  a 
function  of  [02(a)],  at  five  different  pressures  of  N2  in  the  flow  tube. 


[N2]  / 1016  molecule  cm"3 

Figure  5.  Plot  of  the  second-order  rate  coefficient  for  Ca  +  02(a)  as 
a  function  of  N2  concentration.  This  reaction  exhibits  third-order 
(pressure)  dependence  demonstrating  the  formation  of  Ca02( '  A|); 
the  significant  intercept  indicates  that  the  bimolecular  channel  to 
CaO  +  O  is  also  active. 


Figure  4  shows  that 
reaction  1  is  pressure 
dependent.  This  is 
confirmed  in  Figure  5, 
which  is  a  plot  of  k 
against  [N2] .  The  slope 
of  this  plot  yields  the 
third-order  rate 
coefficient  listed  in  Table 
1.  Note  that  there  is  also 
a  significant  intercept  in 
Figure  4,  which  indicates 
that  there  is  a  second- 
order  component  to 
reaction  1. 
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Figure  6  shows  the  first-order  removal  of  Fe  as  a  function  of  [02(a)],  at  two  pressures.  Note  that 
reaction  3  is  clearly  independent  of  pressure.  This  is  confirmed  in  Figure  7,  which  also  shows  the 
second-order  removal  rate  coefficient  for  reaction  2  between  Mg  and  02(a).  The  Mg  reaction  is 
clearly  pressure-dependent  although,  in  contrast  to  reaction  1  (figure  5),  there  is  not  a  significant 
intercept  on  the  ordinate.  The  rate  coefficients  for  reactions  2  and  3  are  also  listed  in  Table  2. 


Figure  6.  Kinetic  plots  showing  the  first-order  removal  rate  of  Fe  as  a 
function  of  [02(a)],  at  two  different  pressures  of  N2  in  the  flow  tube. 
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[N2]  / 1016  molecule  cm"3 


Figure  7.  Plot  of  the  second-order  rate  coefficient  for  Mg  +  02(a)  and  Fe  + 
02(a)  as  a  function  of  [N2].  The  Mg  reaction  exhibits  third-order  kinetics 
forming  Mg02(1Ai).  The  Fe  reaction  shows  no  pressure  dependence,  indicating 
the  formation  of  FeO  +  O  is  the  only  reactive  channel. 
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Table  2.  Rate  coefficients  measured  in  the  present  study  at  296  K.  The  quoted  uncertainties  are  the 
standard  errors  from  kinetics  plots  such  as  Figures  4  and  6.  The  total  uncertainty,  which  mostly  arises 
from  the  systematic  uncertainty  in  the  02(a)  concentration,  is  estimated  to  be  ±40%. 


Reaction 

Bimolecular  Rate  Coefficient 

Termolecular  Rate  Coefficient 

3  11 

cm  molecule'  s' 

cm6  molecule'2  s'1 

Mg  ±  02(a) 

<  2.4  x  10'14 

(1.8±0.2)x  10'30 

Ca  +  02(a) 

(2.7±1.0)  x  10'12 

(2.9±0.2)  x  10‘28 

Fe  +  02(a) 

(1.1±0.1)  x  10'13 

<  2.4  x  10'31 

A  final  point  to  note  here  is  that  Ca,  Mg  and  Fe  react  comparatively  slowly  with  02(X)  [ Campbell 
and  Plane ,  2001;  Helmer  and  Plane,  1994a;  Nien  et  al.,  1993].  Therefore,  even  though  the  ratio  of 
02(a)/02(X)  entering  the  flow  tube  was  in  some  experimental  runs  as  low  as  16%,  the  reactions  of 
these  metal  atoms  with  02(X)  would  have  had  a  negligible  impact  on  their  removal  rates  and  hence 
determination  of  the  02(a)  kinetics. 

Discussion 

The  recombination  reactions  of  Ca  and  Mg  with  02(a)  are  82  and  4020  times  faster,  respectively, 
than  their  reactions  with  ground-state  02(X).  In  order  to  understand  these  striking  differences  in 
reactivity,  we  have  used  electronic  structure  (quantum  chemistry)  calculations,  combined  with  Rice- 
Ramsperger-Kassell-Markus  (RRKM)  theory  (using  the  inverse  Laplace  transform  method  to  solve 
the  Master  equation  [ Campbell  and  Plane,  2001]).  Some  of  this  work  was  done  in  collaboration  with 
Dr  David  Glowacki  of  Bristol  University. 

Figure  8  illustrates  the  stationary  points  on  the  singlet  potential  energy  surface  for  Ca  +  02(a)  (red 
lines),  as  well  as  the  triplet  surface  (black  lines)  which  links  Ca  +  02(X)  with  the  bimolecular 
products  CaO  +  0(3P).  The  hybrid  density  functional  /  Hartree-Fock  B3LYP  method  was  employed 
from  within  the  Gaussian  09  suite  of  programs  [ Frisch  et  al.,  2009],  combined  with  the 
6-311+G(2d,p)  triple  zeta  basis  set..  This  is  a  large,  flexible  basis  set  which  has  both  polarization  and 
diffuse  functions  added  to  the  atoms.  At  this  level  of  theory,  the  expected  uncertainty  in  the 
calculated  reaction  enthalpies  is  ±20  kJ  mol'1.  If  reaction  1  remains  on  the  singlet  surface  then  the 
outcome  is  recombination  to  form  Ca02(*Ai),  a  deep  well  which  is  279  kJ  mol1  below  the  reactants. 
RRKM  theory  models  successfully  the  surprisingly  fast  measured  rate  coefficient  (Table  1)  using  the 
calculated  vibrational  frequencies  and  rotational  constants  for  Ca02,  and  assuming  sensible  values 
for  the  collisional  energy  removal  by  N2.  The  reaction  Ca  +  02(X)  is  much  slower  because  there  is  a 
small  barrier  (~6  kJ  mol'  )  in  the  entrance  channel  and  the  well-depth  of  Ca02(  A2)  is  much 
shallower  [ Campbell  and  Plane,  2001],  as  shown  in  Figure  8. 

The  other  interesting  feature  about  reaction  1  is  the  bimolecular  channel,  which  must  correspond  to 
the  formation  of  CaO  +  O,  which  involves  a  spin  hop  from  the  singlet  to  triplet  surface.  In  an  effort 
to  characterize  the  regions  of  most  likely  spin  hopping  within  the  Ca02  addition  complex,  we  carried 
out  a  number  of  relaxed  multireference  CASSCF  scans  along  the  O-Ca-O  angle.  All  calculations 
were  performed  with  a  10  electron,  7  orbital  active  space.  The  orbital  active  space  consisted  of:  (1) 
the  in-plane  n  and  n*  orbitals  on  O2;  (2)  the  out-of-plane  n  and  n*  orbitals  on  O2;  (3)  the  in-plane  O2 
a  and  a*  orbitals;  and  (4)  the  dz“  orbital  on  Calcium.  The  results  show  that  the  singlet-triplet 
minimum  energy  crossing  point  (MECP)  is  in  the  vicinity  of  the  OCaO  minimum  geometry  (right- 
hand  side  of  Figure  8).  In  order  to  avoid  convergence  problems  during  the  MECP  optimization, 
calculations  on  both  the  singlet  and  triplet  surface  were  carried  out  without  any  point  group 
symmetry.  The  OCaO  singlet  state  was  treated  using  broken  spin  symmetry  in  order  to  capture  its 
diradical  character.  The  results  of  the  MECP  optimization  yielded  a  linear  OCaO  structure  with  Ca- 
O  bond  distances  of  1.83  Angstrom.  Note  that  once  spin-hopping  has  occurred,  there  is  a  large 
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barrier  on  the  triplet  surface  preventing  the  formation  of  Ca  +  C>2(X).  Hence,  this  quenching  channel 
is  closed. 


Ca  +  02  singlet  and  triplet  surfaces 


Figure  8.  Potential  energy  curves  (calculated  at  the  B3LYP/6-311+g(2d,p)  level  of 
theory)  for  the  singlet  surface  (red  lines)  and  triplet  surface  (black  lines). 
Recombination  of  Ca  +  02(a)  produces  mostly  Ca02(1Ai).  However,  there  is  a  conical 
intersection  between  OCaO( 1 A 1 )  and  OCaO(3A2),  where  there  is  a  small  probability  of 
switching  onto  the  triplet  surface  and  generating  the  bimolecular  products  CaO  +  0(  P) 


For  reaction  2,  recombination  to  form  Mg02(1Ai)  is  the  only  possible  product.  The  well  depth  is  149 
kJ  mol1  below  the  reactants,  and  the  measured  rate  coefficient  (Table  1)  is  straight-forwardly 
modelled  using  RRKM  theory.  The  reaction  is  more  than  3  orders  of  magnitude  faster  than  Mg  + 
02(X)  because  there  is  a  significant  barrier  of  23  kJ  mol1  in  the  entrance  channel  of  the  triplet 
surface  [Nien  et  al.,  1993].  Figure  9  is  a  mesh  plot  showing  parts  of  the  singlet  and  triplet  surfaces,  as 
a  function  of  the  distance  between  the  Mg  atom  and  X  (the  mid-point  between  the  O  atoms),  and  the 
Mg-X-0  angle.  There  are  no  crossing-points  between  the  surfaces,  at  any  point  below  the  reactant 
energy  on  the  singlet  surface.  Thus,  formation  of  the  more  stable  MgC>2(  A2),  or  quenching  to  Mg  + 
02(X),  are  ruled  out. 

Reaction  3  between  Fe  and  02(a)  starts  on  a  quintet  surface.  Because  of  the  complex  spin  multiplicity 
of  Fe,  a  detailed  theoretical  treatment  of  this  reaction  involves  surfaces  of  singlet,  triplet,  quintet  and 
septet  multiplicities,  which  is  beyond  the  scope  of  this  study.  Fe(a5D)  +  02( 1  A)  and  FeO(5E)  +  0(3P) 
correlate  on  a  quintet  surface,  explaining  the  observation  of  bimolecular  kinetics  (Table  1).  Initially, 
reaction  3  involves  formation  of  Fe02(5Ai).  This  can  then  rearrange  to  OFeO(5B2),  before 
dissociating  to  FeO  +  O  [Self  and  Plane ,  2003].  However,  Fe02(5Ai)  can  more  easily  dissociate  to 
Fe(a5F)  +  C>2(X).  This  channel  represents  near-resonant  energy  transfer  (A E  =  -10  to  +4  kJ  mol1, 
depending  on  the  Fe  spin-orbit  multiplet),  and  must  be  the  major  product  of  reaction  3.  The  excited 
Fe(a5F)  atoms  will  be  quenched  to  ground-state  Fe(a3D)  in  only  6  ps  at  the  lowest  pressure  of  N2 
used  in  the  flow  tube  [Mitchell  and  Hackett,  1990],  so  this  energy  transfer  reaction  would  not  be 
observable  as  a  loss  of  Fe  atoms. 
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Figure  9.  Potential  energy  surfaces  for  Mg  +  02(a)  (monochrome  shading)  and  Mg  + 
Oi(X)  (coloured  shading),  calculated  at  the  B3LYP/6-311+g(2d,p)  level  of  theory.  The 
diagram  illustrates  that  there  are  no  intersections  between  the  surfaces.  Thus,  the  only 
possible  reaction  of  Mg  with  02(a)  is  recombination  to  Mg02(*Ai). 
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Figure  10  illustrates  the  singlet  and  triplet  potential  energy  surfaces  for  reaction  4,  SiO  +  02(a).  This 
reaction  should  be  extremely  slow,  for  several  reasons:  (1),  a  barrier  in  the  entrance  channel  of  9  kJ 
mol'  ;  (b)  a  barrier  of  55  kJ  mol  in  the  exit  channel  to  Si02  +  0(  P);  (c)  and  the  requirement  to  hop 
from  the  reactant  singlet  electronic  surface  onto  the  product  triplet  surface.  Thus  reaction  4  will  not 
be  competitive  with  the  reaction  SiO  +  O3  [Gomez  Martin  et  al.,  2009b]  as  a  route  for  oxidising  SiO 
to  Si02  in  the  MLT. 
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Reaction  coordinate 

Figure  10.  Potential  energy  curves  (calculated  at  the  B3LYP/6-311+g(2d,p) 
level  of  theory)  for  the  singlet  surface  (red  lines)  and  triplet  surface  (black 
lines).  Note  that  there  is  a  barrier  of  10  kJ  mol1  on  the  singlet  surface  which 
will  slow  recombination  to  OSiCk-  Formation  of  SiCF  +  O  would  require  both 
surface  crossing  and  overcoming  a  significant  barrier  of  55  kJ  mol’1. 


Figure  11.  Vertical  profiles  in  the  mesosphere/lower 
thermosphere  of  the  first-order  removal  rates  of  Si+.  Note 
that  the  major  daytime  removal  process  between  85  and  107 
km  is  by  reaction  with  02(a).  Conditions  are  40°N,  April  at 
midday. 
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Atmospheric  modelling 


The  geophysical  significance  of  these  reactions  was  tested  using  output  from  the  Whole  Atmosphere 
Chemistry  Climate  Model  (WACCM),  which  is  a  state-of-the-art  atmospheric  chemistry  climate 
model  coupled  with  a  detailed  chemistry  scheme  (http://waccm.acd.ucar.edu/).  The  only  reactions 
studied  in  this  project  which  are  atmospherically  important  are  reactions  1  and  5.  Figure  11  shows 
the  vertical  profile  in  the  MLT  of  the  first-order  removal  rates  of  Si+  by  a  number  of  reactions. 
During  daytime,  when  [02(a)]  is  about  30  times  larger  than  [O3],  reaction  5  is  the  most  important 
removal  process  for  Si+.  Between  85  and  105  km,  the  rate  of  reaction  5  is  roughly  double  that  of  the 
reaction  Si+  +  O3  [Eyet  et  al.,  2010].  In  the  case  of  reaction  1,  the  order  is  reversed:  the  rate  of  Ca  + 
O3  [Helmer  et  al.,  1993]  is  roughly  four  times  faster  than  reaction  1,  during  daytime  between  80  and 
110  km. 

Conclusions 

This  was  a  successful  project  in  which  all  major  objectives  were  achieved.  As  explained  above,  it 
was  decided  to  replace  a  study  of  Si  +  02(a)  with  a  study  of  the  Mg  atom  reaction,  both  because  this 
was  experimentally  tractable  (Mg  does  reacts  very  slowly  with  ground  state  02(X),  whereas  the  Si 
reaction  is  close  to  the  collision  frequency),  and  because  the  Mg  reaction  provided  an  illuminating 
contrast  to  Ca  +  02(a).  Although  it  turned  out  that  the  atmospheric  significance  of  these  reactions 
was  limited,  they  have  proved  to  be  of  fundamental  interest  as  a  test  for  advanced  quantum  chemistry 
calculations  and  rate  theory. 
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